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1 Introduction. 

Hard exclusive processes ffl, i.e. involving a large momentum transfer, provide a complimentary 
and an equally important information about the internal structure of hadrons to the one gained 
e.g. in deep inelastic scattering in terms of inclusive parton densities. By means of QCD fac- 
torization theorems physical observables measurable in these reactions, i.e. form factors and 
cross sections, are expressed as convolution of a hard parton rescattering subprocess and non- 
perturbative distribution amplitudes |l| and/or skewed parton distributions || f|, |5|, [| . Moreover, 
it is implied that the main contribution to the latter comes from the lowest two-particle Fock state 
in the hadron wave function. The field-theoretical background for the study of the distribution 
amplitudes is provided by their expression in terms of matrix elements of non-local operators 
sandwiched between a hadron and vacuum states (or hadron states with different momenta in the 
case of skewed parton distributions) 

<K*) = ^Jdz-e**-(0\<p\0)ip{z-)\h). (1) 

Due to the light-like character of the path separating the partons, p, the operator in Eq. ([J) 
diverges in perturbation theory and thus requires renormalization which inevitably introduces a 
momentum scale into the game so that a distribution acquires a logarithmic dependence on it. 
This dependence is governed by the renormalization group which within the present context is 
cast into the form of Efremov-Radyushkin-Brodsky-Lepage (ER-BL) evolution equation |7|, || 

'{x,Q) = V{x,y\a a (Q))®<l>(y,Q), with & = / dy. (2) 



Note that the restoration of the generalized skewed kinematics in perturbative evolution kernel, 
V, is unambiguous and straightforward ||. Therefore, we discuss in what follows only the case 
when the skewedness of the process equals unity. 

Recently, we have addressed the question of calculation of two-loop approximation for the 
exclusive evolution kernels and give our results for the parity-odd sector in Ref . [[Hj . The main tools 
of our analysis were the constraints coming from known pattern of conformal symmetry breaking 
in QCD and supersymmetric relations arisen from super- Yang-Mills theory. In the present note 
we address the flavour singlet parity even case which is responsible for the evolution of the vector 
distribution amplitude. 

2 Anatomy of NLO evolution kernels. 

Our derivation is based on the fairly well established structure of the ER-BL kernel in NLO. Up 
to two-loop order we have 

V(x,y\a s ) = ^V^(x,y) + (£) V^y) + 0(a 3 s ), (3) 

1 



with the purely diagonal LO kernel V^ ' in the basis of Gegenbauer polynomials and NLO one 
having the structure governed by the conformal constraints ITTL W\ O] 



V® = -V ® (V 0) + ^t\ - g®V™ + V (0) ®g + D + G. 



(4) 



Here the first three terms are induced by conformal symmetry breaking counterterms in the MS 

scheme. Contrary to the LO kernel V^ , the so-called dotted kernel V and the g kernel are 

off-diagonal in the space of Gegenbauer moments. They have been obtained by a LO calculation 

T\\ [12|, [1~3|. The remaining two pieces are diagonal and are decomposed into the G kernel 



which is related to the crossed ladder diagram and contains the most complicated structure in 
terms of Spence functions, while the D kernel originates from the remaining graphs and can be 
represented as convolution of simple kernels. 

One of the ingredients of the NLO result are the one-loop kernels. We use for them improved 



expressions of Ref. |I2] which are completely diagonal in physical as well as unphysical spaces of 
moments. In the matrix form we have 



V^ v (x,y) 



Gi 



QQ u{x,y)} -2T F N f QG v v {x,y) 



C F G % v (x,y) 



C A GG v v (x,y)-fS(x-y) 



(5) 



where /3 



§T F N f - ^C A and C A = 3, C F 



4/3, T F = 1/2 for QCD. The structure of the 



kernels reflects the supersymmetry in M = 1 super- Yang-Mills theory |14], [15], [16 



Q v = QQ y a _|_ QQyb 



QG y V ^ QG y a + gQG^c 



GQ y V 



GQy a _|_ 2 G< \ 



GG V V 



<GG„. 



GGyb 



where the functions v l are defined in the following way 



AB„.i 



v\x,y)=6{y-x) AB f\x,y)± 



x — > X 



y->y 



for 



with (here and everywhere x = 1 — x) 



AB f a \ 

ABfbj 
AAfc 

ABrc 



x u(A)-l/2 



y 



v(B)-l/2 



1 

1 
y-x 



x u(A)-l/2 ( 2xy [| _ \ n (xy) 



y 



v(A)-l/2 



+ y — x 



2xy + y — x 



T 



x u(A)-i/2 ( 2x y-x 
y^B)- 1 / 2 { 2xy - y 



for A 



G) 



+ 2 



GG y c 



A = B 



for A 



*B. 



Q 
G 



(6) 



(7) 



The index v(A) coincides with the index of the Gegenbauer polynomials in the corresponding 
channel, i.e. v(Q) = 3/2 and v(G) = 5/2. The eigenvalues of the same f l -kernel in different 



channels are related to each other (here Vjj = Vj 



l QG v a = 6 GQ v a = ^GG y a 



3 6 J j(j + 3) 3 2 ' ( J + i)( J + 2 )' 

,6 = g^6 _ x = _ 2 ^- + 2 ) + 2^(1) + 2, 

QQ V C = _Iqg ! c = 6 Ga,c = Igg f c = 2 

J " 6 ' j(j + 3) j 3 ' J (j + l)(j + 2)(j + 3)- lJ 

Note that we have the identity GC h) c . = QQv°-/3 = GG v®/6, which in the next section will serve as a 
guideline for the construction of V and G kernels. 



3 Construction of V and G kernels. 

To proceed further let us consider first the construction of the so-called dotted kernels whose 
off-diagonal conformal moments are simply expressed in terms of the one-loop anomalous di- 
mensions, ABr ij , of the conformal operators as 6j-2,k{ AB lj — AB lk )djk with djk = — 1[1 + 



-l) jf k ] ( ■ Lr + +k+3 y We introduce the matrix 



Mv, 



C h 



. , ,[ QQ ^(x,y)} + -2T F N f Q G v V (x,y) 
V (x,y) = \ L ,, J + .. , (10) 

! C F GO v V (x,y) C A GG v V (x,y) 



where we use the decomposition analogous to Eqs. (|5]||) for the LO kernels including the same 
"-(-"-prescription although this time the kernels are regular at the point x = y. The general 
structure of AB v l reads 

AB v\x,y) = 6(y-x) AB r(x,y)\n*+A AB r(x,y)±{ X ^*\, for H "j * (11) 

For the dotted a and 6-kernels we have A AB f l (x, y) = with i = a,b. To find the dotted c- kernels 
we make use of the fact that kernels with the same conformal moments in different channels 
are related by differential equations owing to the following simple relations for the Gegenbauer 
polynomials 

&r«* - *> - «#<* - *>■ H§ji c ^* - * - - 6 HtI^(- - >>. <-> 

were w(x\u) = (ii)" -1 '' 2 is the weight function and Nj(u) = 2~ 4u+1 yHvMv+'YX * s ^ ne normalization 
coefficient. From the knowledge of conformal moments, which are determined by the eigenvalues 
of the corresponding kernels given in Eq. (§), and using the expansion of the kernels w.r.t. the 
Gegenbauer polynomials 

ab i, s _ ^ w {x\v{A)) v{A) AB ■ u{B) 

V \ X >y) - 2^ N . ( V (A)) °J+3/2-*/(A)^ L ) V j^j+3/2-u(B)( Z y l ) 



we find then the following differential equations 



GQ-c 



v c (x, y) 



GG-a 



r v a (x,y) + GG v a (x,y) 



d_ 

dy 

4- GQ v\x,y) = -2 QQ v a (x,y)- QQ v a (x,y), 
ax 



QG v c (x, y) 
One of the entries in Eq. (|TTD, namely 



1 d 
3 dx 



GG-c 



'v c (x,y). 



>GGhc 



X 



A GG f c (x,y) = 2-(y-x 



(13) 
(14) 
(15) 

(16) 



has been obtained in Ref. ||13| . Thus defined GG v c kernel possesses the correct conformal moments 



in both un- and physical sectors. Eqs. (|13]jl4]) allow us (after fixing the integration constant) to 
find A G ®f c , while from Eq. (|T5j ) we conclude that /SP G f c is trivially induced by A GG f c . Therefore, 
we have finally 



x 



a gqjc = x 2, 2x _ 3) ]n 1 A QG f c 

y 

Next the g function is given by [|12|, [T 



x 



g(x,y) = 9(y-x) 



f-a 



\ 



in 



(!-f) 



y-x 



C F ^ 



-c, 



3y< 



lnfl-2 

v v 

y—x 



(Ax — 5y + 2xy) . 



(17) 



± 



/ 






y->y 



with (— ) + sign corresponding to (non-) diagonal elements. 

The construction of the diagonal G(x, y) kernel related to the crossed ladder diagrams is 
straightforward up to a number of points which are not obvious and present the most non-trivial 
part of the machinery. Let us give here its construction in more detail as compared to Ref. [|K^ . 



From the result in the flavour non-singlet sector |17|, 0, [19[ and the general limiting procedure to 
the forward case H |§ we know that all entries in the matrix 



G\x,y) 



1 / 2C F (C F 



Qa 

2 



QQG\x,y)\ A 2C A T F N f Q G G\x,y) 
GG G\x,y)] i 



r<2 



2 1 C F C A G QG\x,y) 

must have the following general structure 

AB G\x, y) = 6(y - x) ( AB H l + A AB H l ) (x, y) + 6(y - x) ( AB TT + A AB TT) (x, y), 
with the following expressions for H and H 
AB H\x,y) = 2 



(19) 



(20) 



abA 



AB 



TT(x,y) = 2 



± AB f (Li 2 (x) + Iny Inx) - AB fU 2 (y) 

x\ 1 



(21) 



( AB f T AB T) Li 2 1 



y , 



\n 2 y\+ AB f (Li 2 (y) - Li 2 (x) - lny lnx) 



(22) 



where the upper (lower) sign corresponds to the A = B (A ^ B) channels. For the QQ sector 
we have /S9®H = A^°JI = 0. However, in general these addenda are nonzero and are needed to 
ensure the diagonality of the kernels. From the known two-loop splitting functions we have to 
require as well that in the forward limit these terms contribute only to rational functions and/or 
terms containing single logs of momentum fractions. 

The reduction P — > V B procedure from the forward to non-forward kinematics JT2J is hard 



to handle for the restoration of AH contributions, so we have to rely on different principles. We 

do this by exploring supersymmetry and conformal covariance of M = 1 super- Yang-Mills theory 

15 , |16| . As a matter of fact being wrong for all order results these assumptions hold true within 



the present context since the AB G kernels arise from the crossed ladder diagrams which have no 
UV divergent subgraphs and, therefore, require no renormalization. Thus, these kernels can be 
constructed from six constraints on anomalous dimensions of conformal operators. In principle 
these relations can also be written for the kernels in the ER-BL representation so that taking the 
known two entries of the quark-quark channel one can deduce all other channels. Unfortunately, 
at first glance it seems that not all of these constraints have a simple solution in the ER-BL 
representation. For this reason we modify our construction in the following way. Because of both 
supersymmetry and conformal covariance, the mixed channels are related by the equation 

GQ G\x, y) = ^-Q G GHy, x). (23) 

yy 

Employing this relation we can get a further one, from the so-called Dokshitzer supersymmetry 
constraint, 

^ QQ G l (x,y) + ^ GG G l (x,y) = -3 QG G\x,y), (24) 

which allows to obtain the GG entry provided we already know the kernel in the mixed channel. 
Let us consider first the parity odd sector. At LO we have for moments Q G v A = 6®%? Thus 
we can simply obtain the QG kernel differentiating the QQ one. Fortunately, it turns out that 
the Q G G A kernel can be obtained^] in the same way 

Q G G A (x,y) = ^G a (x,y), (25) 

dy 

where QQ G a is given by Eqs. (H-H) with A QQ H a = A QQ H a = 0. The GQ entry is simply deduced 
from Eq. (|23"D , while the GG one comes from the solution of the differential equation ([24]). The 



integration constant as a function of y is almost fixed by the necessary condition of diagonality 

GG G A (x,y) = ^ GG G A (y,x). (26) 

2 The correctness of this and subsequent results is checked by forming the Gegenbauer moments and comparing 
them with known NLO forward anomalous dimensions [BCI . 



The remaining degree of freedom can be easily fixed from the requirement that the moments 
GG Gj i are diagonal for i = 0,1. To simplify the result, we remove a symmetric function (w.r.t. the 
simultaneous interchange x — > x and y — > y) which enters in both A GG H and A GG H kernels, 
however with different overall signs and, therefore, disappears from GG G . We present our final 
results in a symmetric manner as (cf. [pi 



A QQ H A (x,y) = A QQ H A (x,y) = 0, (27) 

A^ G H A {x,y) = A^ G H A (x,y), A^ G H A (x,y) = -^-A G ^H A (y,x) (28) 

{yy) 

A GQ H A (x,y) = -A GQ H A (x,y) A GQ H A (x,y) = -2— Ins + 2— Iny, (29) 

y y 

A GG H A (x,y) = -A GG H A {x,y), (30) 



1 — 2xx 1 — 2a; (1 + x) x(x + y — 3xy) x(x + y — 3xy) 

4^2 ^2 yy2 yy2 



Now instead of dealing with the whole parity even sector, we can consider only the difference 
between vector and axial-vector functions 



H v =H A + H d . (31) 

In LO we have the simple equation GC hj c . = Q°v^/3 = GG t>"/6, see Eq. (£5|), which allows us to write 
down a simple relation between the kernels in different channels. However, to preserve the generic 
form of the GO -G c function in the forward limit [20| we have used the following modified differential 



equations^: 



^_GQ H c _ _ 4 (qQjjo, _j_ g QQrc ^ QQ /-A ^_GQ-^c _ 2 /GGjya + g GG/^c ^ GG« 

dx V J ' dy \ 



'' n w = -4 ( QQ r + 9 QQ r I QQ r) , — gq h c = 2 ( GG r + 2 GG r l GG r) , (32) 



dy 

d_ 

dx V ) ' ' dy 

where f c (x,y) = f c (x,y). The kernels GG H a and GG fT are the parts of the whole parity odd 
functions derived in the fashion already explained above. The two sets of differential equations 
can be solved up to two integration constants which can easily be fixed from the diagonality of their 
conformal moments. Finally, we simplify the solution by adding pure diagonal pieces containing 
a and c kernels and their convolution as well as by removing symmetric terms which die out in 
G ®G. 

The entry in the QG channel can be obtained from the supersymmetric relation ( p3|) . To 
construct the GG kernel we use then the constraint ( |2J|) with °-QG c = 0. We determine the 



integration constant as a function of x in the same manner as described previously. Our findings 



3 Note that we introduce a shorthand notation for the convolution, namely, ®®f l <g>° QQfi is understood as 
convolution of the corresponding ER-BL kernels and then keeping only the part proportional to 9(y — x). 



r r 

for the A AB H and A AB H can be summarized in the formulae 



A QG H\x, y) 



(yy) 2 

A GQ H\x,y) = A GQ H\x,y) + 20 



A G ^H\y, x), A^ix, y) = ^-A^y, 



iyyf 



-arrrS,- ■, ~„x(% _ %) ^(3 + 2x) _ x(3 + 2x) 



3y 



32/ 



lnx + 4- 



3y 



Iny, 



A GQ H 5 {x,y) 



61 



+ 2xx(l - (3 - lOx) lny + (3 - lOx) lnx 



x (6 — 19x + 6x 2 ) „x (y + x(x — x)) , x (y + x(x — x)) , 
— i-2^^ ^ — lny + 2^^ —Ins, 



32/ 



GGWf- 



A GG tf (x,2/) = A GG H (x,y) 



y y 

20 - 18x + 55xx 20 - 23xx 17 + 32x + 28x 2 



6y 2 



6y 2 



62/2/ 



xlx — x 2 + 3a; \ x ( x — x 2 + 3s\ 
- 2- [ 2^— + Jinx -2- (2 + — — ] my, 



n i/ 



2/ 



n y 



</ 



A GG Il\x,y) 



, /20-22x + 21xx 20-22x + 21xx 39 + 38xx" 
■(1 - x - y) — + — + 



6r 



63/ 



•i iZ 



62/2/ 



/x 3 x 2 (21-20x) xx 2 \ / x 3 x 2 (21 - 20x) xx 2 ' 
+ 2 — K —- L -2— lnx + 2 — ^- 1-1- 



(33) 
(34) 



(35) 



M 



32/ 



y 



3y 



V 



In 2/. 



4 Restoration of remaining diagonal terms. 



As in the twist-two axial and transversity sectors |T(J it turns out that the remaining diagonal 
piece, D , can be represented as the convolution of simple diagonal kernels. To find it we take 
first the forward limit 



1 / QQy IQGy 
P(z) = LIM V(x, y) = lim — 

V 7 V ' t-*Q\t\ \ zGQy IGGy 



z 1 



T T / 



and compare our result with the known DGLAP kernel P [20|. In this way, 

D v (z) = P v (z) - LIM l-V I (v^ v + ^t) - g § V^ v + V (0)y ® 5 + G l 



(36) 



we extract the remaining part LIMZ) y (x, 2/) and find then the desired convolutions in the forward 
representation. Note that we map the antiparticle contribution, i.e. z < 0, into the region z > 
by taking into account the underlying symmetry of the singlet parton distributions. Our findings 
can be immediately mapped back into the ER-BL representation: 



QQjf 



QGjf 



Cl [D F ], - cA [D ], - C F ( <?,- 



C 



\Q% + 2 Q % a 



+ AC F T F N f {- 
= -C F T F N f \2 



a « nn a 



QQ v a _ 6 Q V I QQ V C _ Q% a + Z$V 

6 

(§) Q( V + Q< V ® Q< V + - Q( V + 6 Q< V 



(37) 



(38) 



GQ D ] 



GGjf 



+ 2C A T F N 



f 



QQ, 



+ 56 Q( V 



,QG, 



130 



a e rtH o. 



cl 



— Cpf3 
+ CfCa 



+ 



GG v 



GG y 



? " 2C < 2 > 



® 



:«Q,, a 



QG« 



GQ, 



QG V ° 



'■v + 3 G S 



GQ„. _i_ CQ, 

2 



' v + 



-5 G V® G V-3 G V 

+ - GG y (g) Gg u + -H; 
4 3 



GG / 



G 



V 



- G V 



6 



, G V 



43 

9~ 



- + 2C(2)) G V 



^ 



GG V A 



2 

+ 3 

- C 



GG y A 



GG a , HgG v c 

3 



131 gg° , !)i ^v ' 



+ g-4 C (2)) G v}, 
-14 G V® G V + 12 G V® G S C 



/3o 



. +— GG v -2S(x-y) 

12 18 v j; 

5 



._GG (g,^ _|_ _ 

2 3 



G S V 



13 



GG„. a 



f3 ^V + ^ G V + 25(x-|/) 



+ CVTVA^ j G V ® G V + ^ G V -5(x-y)\. 



(39) 



(40) 



where Dp, Dp functions are known from the flavour non-singlet case (|IT]]. In comparison to the 
parity odd sector the convolution of c-kernels appears as a new entry. It is worth mentioning that 
our result for the evolution kernels in the parity even singlet sector possesses the correct conformal 
moments in both the physical and unphysical sectors. This is to be contrasted with an explicit 
momentum fraction space calculation at LO and quark bubble insertions in NLO kernels for the 
mixed channels [12| where the improved kernels do not appear. 



5 Conclusions. 



To recapitulate the results presented here, we have reconstructed the two-loop singlet evolution 
kernels responsible for the scale dependence of the vector meson distribution amplitudes. We 
have avoided cumbersome next-to-leading calculations by adhering to an extensive use of the 
conformal and supersymmetric constraints derived earlier which thus play a paramount role in 
the formalism. The correctness of the results given here is proved by evaluating the Gegenbauer 



moments of the kernels which coincide with the anomalous dimensions derived in Ref . |13| . Our 



findings allow to use now the direct numerical integration (see Ref. [21, 22 1 for a leading order 
analysis of non-forward parton distributions) of the generalized exclusive evolution equations which 
provides a competitive alternative to the previously developed methods of orthogonal polynomial 



reconstruction of skewed parton distribution pursued by us in Ref. |23| . 
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